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Abstract 



The direct charmonium and bottomonium production rate observed at high- 
p± in pp collisions at the Tevatron is factors of ten larger than predictions 
based on conventional perturbative QCD. We show that this excess can be 
accounted for by our model for soft colour interactions, previously introduced 
to describe in a novel way the large rapidity gap events observed at HERA. 
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Heavy quarkonia, i.e. bound states of heavy quark-antiquark pairs, are thought to pro- 
vide a useful testing ground for perturbative Quantum Chromo Dynamics (pQCD). The 
reason being that the charm and bottom quark masses provide a large scale which makes 
the short distance behaviour calculable using a perturbative approach whereas the nonper- 
turbative contributions can be factorised into a wavefunction. It was therefore a surprise 
that measurements by the CDF and D0 collaborations at the Fermilab Tevatron pp 
collider (a/s = 1.8 TeV) gave cross-sections of direct high-p^ J/ip^ip' and T, T' production 
far above the expectation from the colour singlet model based on conventional pQCD. 
The observed excess is generally an order of magnitude and increases to a factor ~ 50 with 
increasing p^. 

In the colour singlet model pQCD is used to calculate the production of a QQ pair in 
a colour singlet state. This forms a quarkonium state with the same angular momentum 
quantum numbers "^^^^Lj by coupling to the non-perturbative wave function at the origin, 
which is obtained from the leptonic decay width. The striking failure of this model has 
led to several phenomenological investigations and some new models; for a review see 0. 
The colour octet model, which is based on non-relativistic QCD, also takes into account the 
more abundant perturbative production of QQ pairs in colour octet states. The unknown 
probability for the transformation into colour singlets due to non-perturbative processes is 
parametrised with matrix-elements that have to be fitted to data, but are universal so that 
the model can be tested by studying several different processes. Similarly, in the colour 
evaporation model a certain fraction of all QQ pairs, independently of their production 
process and quantum numbers, form a quarkonium state. Thus, both the colour octet model 
and the colour evaporation model requires fitting to experimental data. In contrast, our soft 
colour interaction (SCI) model gives a prediction also of the absolute rate which, as will 
be shown in this Letter, is in good agreement with the Tevatron data. 

The SCI model was introduced as a novel way to interprete the rapidity gap events 
observed in deep inelastic scattering (DIS) at HERA. The model does indeed describe the 
salient features of the ZEUS and HI data. The conventional interpretation of these 
events is in terms of deep inelastic diffractive hard scattering ^ on partons in the pomeron, 
which is a colourless object exchanged from the quasi-elastically scattered proton. Although 
phenomenological models based on this idea work quite well to describe the data, there 
are theoretical problems with the pomeron approach. In particular, the factorisation into a 
pomeron flux, a pomeron parton density and a hard interaction described by a QCD matrix 
element may not be universal for all processes, e.g. DIS and hadron collisions. The pomeron- 
proton interaction is soft and thereby occurs on a long space-time scale. It may therefore be 
incorrect to consider the pomeron as decoupled from the proton during and after the hard 
scattering. 

It is more natural to expect soft interactions with the proton both before and after the 
snapshot of the DIS probe. To investigate this line of thinking, we have developed a model 
with a mechanism for soft colour interactions as an alternative to the approach based on 
the pomeron and Regge phenomenology. 

The basic idea in our SCI model is that there may be additional soft interactions, not 
previously accounted for, at a scale below the cut-off Qq for perturbative QCD. Obviously, 
interactions do not disappear below this cut-off. The question is how to take them into 
account properly and whether conventional hadronization models give a complete descrip- 
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FIG. 1. Main pQCD processes for production of a QQ pair in (a) leading order and (b) 
next-to-leading order. 



tion. We propose that the quarks and gluons produced by conventional pQCD processes, 
as described by matrix elements and parton showers, interact non-perturbatively with the 
colour medium of the proton remnant. This should be a natural part of the process in which 
'bare' perturbative partons are 'dressed' into non-perturbative ones and the formation of 
the confining colour flux tube in between them. These soft interactions cannot change the 
momenta of the partons significantly, but they change their colour and thereby affect the 
colour structure of the event. This in turn will lead to a modified hadronic final state, e.g. 
with rapidity gaps or quarkonium. 

Lacking a proper understanding of non-perturbative QCD processes, we have constructed 
a model P| to describe and simulate such soft colour interactions within conventional Monte 
Carlo (MC) event generators. All partons from the hard interaction plus the remaining par- 
tons in the proton remnant are assigned appropriate colour and anticolour charges. Partons 
from the hard interaction can make a soft interaction with partons in the proton remnant 
exchanging their colours but not changing their momenta, which may be viewed as soft non- 
perturbative gluon exchange. Similarly, partons in the proton remnant can exchange colour 
whereas no direct colour exchanges are allowed between the perturbative partons. Since the 
process is non-perturbative the exchange probability for a pair cannot be calculated from 
first principles and is instead described by a phenomenological parameter R. The number of 
soft exchanges will vary event-by-event and change the colour topology. In the Lund string 
model [^] this corresponds to a modified string stretching that alters the outcome of the 
hadronization. 



The SCI model was first implemented in Lepto |[T1]] simulating DIS. Here, the main 
effect producing rapidity gap events at HERA is SCI between the partons from the hard 
scattering and those in the proton remnant. In particular, with gluon initiated processes at 
small Bjorken-x, the colour octet charge of the produced hard parton system {qq plus possible 
additionally emitted partons) can turn into a colour singlet. This may give a hadronic system 
X separated in rapidity from a very forward moving proton remnant system of small mass 
(e.g. a proton or a resonance). The main features of diffractive scattering emerges naturally 
and the HERA rapidity gap events can be surprisingly well described by this simple model 
i- 

This SCI model can be straightforwardly applied to describe heavy quarkonium produc- 
tion in pp] a colour octet QQ pair from pQCD can be turned into a singlet and thereby form 
an onium state, provided that the mass is appropriate. In leading order (LO) pQCD heavy 
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quark production occurs through gg — > QQ (Fig. |T|a) and qq QQ. However, higher order 
processes involving gluon sphtting g QQ are important. For example, the next-to-leading 
order (NLO) process gg gQQ illustrated in Fig. |I]b gives a large contribution because it is 
an as correction to the large cross-section for gluon scattering {gg gg). Matrix elements 
with explicit heavy quark masses are available up to NLO. However, since the virtuality 
of the gluon need not be very large to split, in particular to cc, still higher orders may be 
important at Tevatron energies. These can be taken into account approximately through 
the parton shower approach. 

The LO and parton shower production of QQ pairs are available in the MC generator 



Pythia [jT2[. On the generated parton level events we apply the above SCI mechanism 
which will turn some QQ pairs into colour singlets. A quarkonium state is then produced if 
the invariant mass rnqq is below threshold for open heavy flavour production (2m a/). Thus, 
the cross-section is 

r2mM da^ 

CTonium = ^"^QQ (l) 

J2mQ amQQ 

where the singlet QQ cross-section da^ = da ® SCI is obtained from the application of the 
SCI model on Pythia events. This is not just a constant fraction of the original QQ cross- 
section, but depends somewhat on the partonic state and the possible string configurations. 
Whereas the number of possible string configurations increases with parton multiplicity, the 
number of string configurations giving a singlet QQ pair is more or less constant. Therefore, 
the relative probability for a singlet QQ decreases slowly with parton multiplicity. This 
in turn, implies a slight decrease (~ 10 % in the observed p^-range) of the singlet QQ 
fraction with increasing in the hard scattering process which causes more abundant 
parton showering. 

The invariant mass of the QQ pair is in principle given by the pQCD process. How- 
ever, we expect this mass to be smeared by the soft colour interactions involving energy- 
momentum transfers of the order a few hundred MeV. Irrespective of the original QQ mass 
(below 2m jv/), we therefore divide the quarkonium cross-section onto the different quarko- 
nium states based simply on spin-statistics as suggested by [Q. This should be a good 
approximation since the heavy quark system is nearly non-relativistic. The cross-section for 
a given quarkonium state X with total angular momentum Jx is then given by 

^X ^ onium (2) 

where Vx = (2Jx + ^)/nx corresponds to a partial width. Here, we have included a sup- 
pression of radially excited states, i.e. with the main quantum number nx- This reproduces 
approximately the differences between different if) and T states regarding their leptonic width 
and thus the wave function at the origin |T^, which should be particularly relevant when 



the essentially pointlike QQ pair forms an onium state. 

The results of this model are shown in Fig. |^ and Table |I| together with the Tevatron 
data. The agreement is remarkably good, considering the simplicity of the SCI model and 
the fact that it was originally contructed for a different physics issue. As mentioned, the 
absolute normalization is not adjusted to data, but given by the model. The only parameter 
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in the SCI model, i.e. the probabihty i? of a colour exchange between two partons, is kept 
at the value R = 0.5 chosen to reproduce the rate of rapidity gap events in DIS at HERA. 
However, as for the rapidity gap rate, the dependence of the onium rate on this parameter is 
quite small. For example, changing to i? = 0.1 only decreases the cross-section with ~ 30 %. 

There is however a strong dependence on the heavy quark mass mq. We have used 
rric = 1.35 GeV and rrib = 4.8 GeV. Increasing the quark masses by 0.25 GeV decreases the 
cross-section with a factor 3, whereas a decrease of the quark masses by the same amount 
increases the cross-section with a factor 2. This is as usual in pQCD heavy quark production 
close to threshold and can be understood from Eq. (p. 

The cross-section also depends somewhat on other features of the Monte Carlo model like 
the QCD A parameter, details in the parton shower model and the choice of parton density 
parametrizations. These issues are not particular for the SCI model, but of a general charac- 
ter. We therefore do not adjust any of these features, but leave them as they are by default 
in the Pythia MC model which provides a standard for many processes and observables. 
We note that there are no significant changes in our results when varying the virtuality scale 
for the parton shower or using different available parton density parametrizations (e.g. MRS 
or CTEQ) based on recent data. However, including the multiple parton-parton interactions 
in Pythia lowers the total onium cross-section slightly (~ 20 %). 

The observed p± spectra of the different onium states are also quite well reproduced by the 
model, as demonstrated in Fig. |^. The contribution from the LO process (shown separately 
as dashed curves) have a steeper p± dependence and only contributes significantly at lower 
p±, in particular for Thus, the higher order contributions are most important. 

Comparing model and data in detail, there is a tendency for the model to have a slightly 
flatter p± distribution for the J/ ip production. This may be an effect of using the approxi- 
mate parton shower gluon splitting g ^ cc instead of the NLO matrix element which might 
be more appropriate at large p±. It may also be influenced by the detailed contribution of 
different onium states, since J/ip from decays contribute more at smaller p±. 

The relative rates of J/ from direct production and from decays of higher onium states 
are given in Table | together with some ratios of cross-sections for different onium states. 
The model gives good agreement with data in most cases. This shows that simple spin 
statistics can describe the main effects. The only deviation of some significance is the ratio 
cr(T(2S))/(T(T(lS)), where the simple suppression of radially excited states seems too strong. 

The SCI model provides a mechanism for how a colour octet QQ pair is turned into 
a singlet and thereby an absolute normalization of the cross-section for the quarkonium 
production. In addition, one obtains a direct relation to other manifestations of soft inter- 
actions, such as diffraction and rapidity gaps. In the colour octet model the probability to 
form a singlet is parametrized into matrix elements, but no explicit mechanism for how an 
octet is turned into a singlet is given. 

These models are in clear contrast to the colour singlet model for quarkonium production 
which cannot explain the large rate observed at the Tevatron. In addition, they differ in 
the fragmentation function of J/ip, i.e. Dg^j/^{z, fi) where z is the momentum fraction of 
the J/ip relative to the 'mother' gluon and /i is the factorisation scale. In the colour singlet 
model the fragmentation function is relatively fiat whereas in the colour octet model it 
is essentially a (5-function at ^ = 1 |jl5| . 



The fragmentation function can be measured by taking the ratio of the transverse mo- 
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FIG. 2. Transverse momentum spectra of prompt charmonium J/tp, ip' and bottomonium T, T' 
in pp collisions at the Tevatron energy. Data from CDF [|| (filled symbols) and D0 Q (open sym- 
bols) compared to results from the soft colour interaction (SCI) model applied to cc/bb production 
from leading order (a^) matrix elements (dashed lines) and with the inclusion of higher order 
contributions calculated in the parton shower approach (full lines). Production central in pseu- 
dorapidity (a) |r/| < 0.6, (b) \y\ < 0.4 and forward (c) 2.5 < \t]\ < 3.7 (where an assumed 20 % 
systematic error has been added in quadrature). 
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TABLE I. Relative rates in % for different onium states in the Monte Carlo model and the 
Tevatron data Q. The charmonium cross-sections are in general for p± > 5 GeV and \r]\ < 0.6 
whereas the bottomonium cross-sections are for \y\ < 0.4. 

MC model Data 

J/^p direct 61 64 ± 6 

J/ip from Xc 26 30 ± 6 t 

J/tp from tp' 13 7 ± 2 ^ 15 ± 5 t 

^(X^)/(^(XJ) + ^(X^)) 62 47 ±9 

a{ip')/a{J/i:) 32 25 ± 6 § 

cj(T(2S))/a(T(lS)) 29 46 ± 9 

^ Data are for p_L > 4 GeV. 

^ Lower (upper) value for ~ 5 (18) GeV. 

^ Assuming the same fraction of J/^ and ip' from h hadron decays. 



mentum of the J /ip and the total transverse momentum in a cone centered at the J /ip ||16|| . 
Following conventional jet algorithms for high energy hadronic collisions, one may choose a 
cone Ai? = a/A0^ + Ar]'^ ~ 0.7 in pseudorapidity and azimuthal angle. Applying this to our 
SCI model results in a fragmentation function for J/ip's with p± > 5 GeV which is strongly 
peaked at ;z = 1. Although its exact form depends on details in the model, it is essentially 
as expected in the colour octet model and very different from the colour singlet model. 

In summary, we have shown that the large rate of direct high-p^ charmonium and bot- 
tomonium in high energy pp collisions can be quite well described the Soft Colour Interaction 
model. The same SCI model also accounts for the rapidity gap events observed at HERA. 
This indicates that these features of non-perturbative strong interactions are quite general 
and can be described by a universal model. 
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